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We present the experimental observation of a magnetically tuned resonance phenomenon resulting 
from the spin mixing dynamics of ultracold atomic gases. In particular we study the magnetic field 
dependence of spin conversion in F=2 *^Rb spinor condensates in the crossover from interaction 
dominated to quadratic Zeeman dominated dynamics. We discuss the observed phenomenon in the 
framework of spin dynamics as well as matter wave four wave mixing. Furthermore we show that 
the validity range of the single mode approximation for spin dynamics is significantly extended in 
the regime of high magnetic field. 
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Ultracold quantum gas spin mixtures are currently re- 
ceiving rapidly growing experimental and theoretical in- 
terest. They combine the unprecedented control achieved 
in single component Bose-Einstein condensates with in- 
trinsic degrees of freedom. The spin-dependent coupling 
connects quantum gas physics to fundamental magnetic 
phenomena. Earlier work analyzed the basic magnetic 
properties of these spinor condensates P, S S 01 and 
demonstrated antiferromagnctic behavior in the F = 1 
state of 23Na ^ and the F = 2 state of ^'^Rb [i| as well as 
ferromagnetic behavior in the F = 1 state of *^Rb Q ■ 

A particularly interesting feature of these clean and 
undisturbed systems is that they give access to quantum 
aspects of magnetism such as interaction-driven spin os- 
cillations p, H, 1^ 0, 0, 0, , the quantum classical 
transition [13 or domain formation 13 ■ The mag- 

netic field dependence of the underlying spin conversion 
process has been analyzed in detail for F = 1 spinor con- 
densates [ij,[ia,ll3ll3 and besides oscillatory behavior a 
resonance phenomenon was predicted ^3- Experiments 
so far focused on magnetic fields above the resonance ^3 
and on magnetic phase controlling of the spinor dynam- 
ics |l2j but a spin mixing resonance has not yet been 
observed. 

In this letter we study for the first time magnetic field 
tuned spin mixing in _F = 2 ^^Rb spinor condensates. 
In particular as a main result of this letter we demon- 
strate resonant spin mixing behavior in spinor conden- 
sates. This resonance phenomenon becomes evident in 
the temporal population evolution shown in Fig. ^ As 
a surprising central effect, the mean field driven spin dy- 
namics (population transfer from mp ± 1 to mp — 0) 
strongly depends on the external magnetic field. In con- 
trast to studies on the natural ground state phase of 
the system and early findings on spin dynamics suppres- 
sion 11,0 at high magnetic field values we now find, that 
an external magnetic field not only leads to adverse ef- 
fects but can even stimulate spin dynamics. The spin 
mixing amplitude is not large for zero magnetic offset 
field, but shows a pronounced maximum for finite mag- 
netic field values as predicted in The observation 
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FIG. 1: Left: Temporal evolution of the mp = 0-population 
for different ratios of quadratic Zeeman and interaction energy 
(see text). Right: Magnetic field dependence of the amplitude 
and period of the mF = population as extracted from a 
sinusoidal fit to the oscillatory part of the data. 



and understanding of this resonance in self driven spin 
dynamics in quantum gases is central for detailed mag- 
netic state manipulation. A q uantum control for general 
spin systems (for F=l see [ij) based on this spinor res- 
onance process might lead to new schemes for spin en- 
tanglement and quantum information applications with 
spinor gases. 

The experimental setup has been described in |^ and 
the procedure resembles the one in We prepare de- 
generate ^^Rb ensembles in the F = 2, nip = —2 state 
containing typically 3 • 10^ atoms in an optical dipole 
trap (trapping frequencies uJx ■ i-^y '. uj^ ~ 2tt ■ 14, 5 : 
105 : 600 s~^). A radio frequency 7r/2-pulsc of typi- 
cally 40/iS duration rotates this state to our initial state 
C™(0) = (1/4,1/2,^6/4,1/2,1/4)^ I23. After hold- 
ing the sample for a variable time, t, at a well defined 
magnetic offset field we switch off the trapping potential. 
The population in each spin component is extracted from 
absorption images after Stern-Gerlach separation during 
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time of flight |^ . Fig. shows the temporal evolution 
of the mp = 0-population in dependence of the relative 
value of the quadratic Zeeman frequency, q and the mean 
field spin mixing frequency, gi (n) . The magnetic field de- 
pendence is parameterized to first order hy q = p^/whfs 
with p = hbB /2h and the *^Rb ground state hyperfine 
splitting ojms ~ 27r- 6.835 GHz. The spin dependent mean 
field interaction is proportional to the density, ri, of the 
sample and for ®^Rb dominated by the parameter gi in- 
troduced below. 

Our observations show that the whole spin mixing dy- 
namics depends critically on the relative size of the en- 
ergies associated to the quadratic Zeeman shift and the 
spin dependent mean field interaction. Interestingly the 
oscillation amplitude in Fig. ^ is reduced whenever one 
of these energies dominates, i.e. in both the quadratic 
Zeeman regime {q » gi{n)) and the mean field regime 
{q « gi{n)). A maximum occurs close to balanced 
quadratic Zeeman and mean field contributions. The im- 
portant point is, that these two effects compete with each 
other. The quadratic Zeeman effect lifts the energy of the 
mp = state with respect to the other states, while the 
mean field interaction lowers it. More precisely the phase 
evolution due to these energy contributions is occurring 
in opposite directions. However, as spin mixing is a coher- 
ent process, this phase determines the direction of spin 
mixing. The accumulation of many particles in one state 
can thus only happen for slow phase evolution, i.e. close 
to the balanced case. 

We want to emphasize, that the many-body spin mix- 
ing resonance presented in this letter crucially relies 
on nonlinear many-particle effects like self- and cross- 
phase modulation. It is therefore fundamentally differ- 
ent from the " Rabi" -resonance case observed in spin- 
changing two-particle collisions i^i which these ef- 
fects are absent. 

In the following we will begin our analysis of the ob- 
served spin mixing resonance phenomenon with a discus- 
sion of the the limiting cases. In particular we present 
and interpret analytic solutions for the asymptotic be- 
havior in these regimes. Our analysis is based on the 
evolution due to the mean field hamiltonian developed 
in 3 using the single mode approximation (SMA). 

Spin mixing occurs as a consequence of two-particle 
collisions and can be classified by the total spin / of 
the colliding pair. For F = 2 the allowed values are 
/ = 0,2,4. The spin mixing interactions are parame- 
terized by the coupling coefficients gi = "4-02 g^j^j 
32 = i^| 7ao-i0a2+3a4 [| ^ ^gj^^g ^^^^ spin-dcpendcnt 
s-wave scattering lengths a/. The dynamics of F = 2 
spinor condensates will in general be a complex super- 
position of the different mixing processes summarized in 
table CK we use the notation (i) for waves/particles in the 
mp = i state). This table also lists the relative quadratic 
Zeeman energy shifts of the involved states. ®'^Rb has the 
advantage that the value of 32 is very small , such that 
the processes listed in the lower part of table HI is largely 
suppressed and can be neglected in the following [23| ■ In 



TABLE I: Four wave mixing processes in F = 2 spinor BEC. 



process 


coupling 


q.Z. energy difT. 


(0) + (0) ^ (+1) + (-1) 


yi 


2a 


(+1) + (+1) ^ (+2) + (0) 


91 


4q 


(-1) + (-1)^(0) + (-2) 


91 


4q 


(+l) + (-l)^(-f2) + (-2) 


91 


6q 


(0) + (0)-(+l) + (-l) 


92 


2q 


(+l) + (~l)^(+2) + (-2) 


92 


6q 


(0) + (0) ^ (+2) + (-2) 


92 


8q 



addition we will concentrate on the specific initial state 
C'"'(0) used in our experiments. For the population evo- 
lution in the limit of low magnetic field (<? << gi{n)), i.e. 
the "mean field region", we find to first order: 



^ ' 1 - -^-^ \cos{Agi{n)t) - 1 
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2,91 (n) 
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2.91 (n) 



[cos(45i(n)0-l] 
- [cos{4gi{n)t) ^ 1] 



(1) 



This mean field dominated evolution is just a sinusoidal 
oscillation between the spin components and shows a 



magnetic field independent period T^f = '^gi{n) 
served in Fig. ^ Furthermore these equations indicate 
that the oscillation amplitude in this regime should be 
proportional to the quadratic Zeeman shift g, which is 
reflected in Fig. ^ 

In the other limit at high magnetic field, i.e. the 
quadratic Zeeman regime with q >> gi{n), the popu- 
lation evolution is dominated by the quadratic Zeeman 
energy and we find: 
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The oscillation amplitude in this regime decays with i 
in accordance to the data in Fig. ^ However now sev- 
eral frequencies are involved as a result of the different 
mixing processes in F = 2 spin dynamics. The os- 
cillation periods refiect the different quadratic Zeeman 
shifts for the possible gi mixing processes as listed in 
tabled For the chosen initial state, spin mixing is dom- 
inated by the cos{2qt) term in Eq. |21 connected to the 
(0) -I- (0) ^ (—1) -I- (-1-1) process, which has a factor of 
four higher weight than the cos(4g<) term. 

From the above tendencies in the population oscilla- 
tion, it is clear, that a maximum amplitude must occur 
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at an intermediate magnetic field value, i.e. the reso- 
nance found in Fig. ^ 

In the following we will develop some more qualita- 
tive physical insights by an interpretation of the reso- 
nance phenomenon in terms of " phase matching" in four 
wave mixing (see Fig. First we note that given the 
prefactor the oscillation amplitude of the mp ~ ±2 
wave is relatively small in both the mean field regime 
and the quadratic Zeeman regime. This justifies an 
approximate view as a pure (—1) -I- (+1) ^ (0) -I- (0) 
four wave mixing process. In contrast to nonlinear op- 
tics or four wave mixing in single component conden- 
sates ^3 we do not deal with wavevector or momen- 
tum modes |23, |^ but with spin modes, i.e. here the 
mp = -fl wave and the mp = wave couple to two 
times the mp = Q wave (this view is justified, as we 
are considering trapped samples in a single momentum 
state). Spin mixing in trapped samples corresponds to 
degenerate collinear four wave mixing |2ll | and (in con- 
trast to single component four wave mixing) in principle 
allows infinite interaction times, only limited by finite 
temperature effects H [Illll iH Hi ■ 

In view of the optical analogy it is obvious, that phase 
matching considerations are essential to understand the 
resonance in spinor four wave mixing. The important 
point is, that the value of the relative phase of the ini- 
tial and final waves/spin components determines the di- 
rection of wave mixing. In our case the (0) wave will 
be populated if the combined phase of the (-1-1) + ( — 1) 
waves is ahead of twice the (0) wave phase, i.e. 9 = 
(f+i + ip-i — 2(^0 G [0..7r] and it will be depopulated, if 
6 G [k..2it] (modulo 27r). 

In general the relative phase evolution of the spinor 
components in spin mixing is highly nontrivial, as in ad- 
dition to the quadratic Zeeman energy shifts it depends 
on the spin coupling and the spin component populations 
in a nonlinear way. We find that for our system, the com- 
petition between mean field driven dephasing (tending to 
decrease 9) and quadratic Zeeman shift driven dephasing 
(tending to increase 9) determines the evolution of the 
system. Most importantly the evolution depends on the 
relative size of the quadratic Zeeman energy shift and 
the maximally achievable mean field shifts, as shown in 
Fig. 121 For large magnetic fields, i.e. always negligi- 
ble mean field energy shifts, 9 will continuously grow, 
i.e. the (±1) wave evolves faster than the (0) wave, and 
the population transfer shows an oscillatory behavior de- 
pending on which wave is lagging behind at which instant 
(Fig.|2Ji). In this regime the oscillation period is expected 
from Eq.^to be given by ^. Furthermore the oscillation 
amplitude, which depends on the unidirectional mixing 
time, should be proportional to this period. Both these 
expectations are confirmed by the data, as shown in the 
inset of Fig. 

For small magnetic fields the mean field energy shift 
grows with increasing population in the (0) wave, until 
it exceeds the quadratic Zeeman shift and thus reverses 
the evolution of 9, i.e. the (0) wave will speed up and 




FIG. 2; Schematic representation of four wave mixing in a) 
the quadratic Zeeman regime and b) the mean field regime. 
The direction of spin mixing is indicated in the upper and 
lower half of the circle symbolizing the relative phase 6. For 
6 > Q the (0) waves (upper waves) get populated, while for 
6 <Q the (+1) and (—1) wave (lower waves) increase. 



decrease its lag behind the (±1) wave, eventually overtak- 
ing it. In this case 9 will remain confined in the interval 
[— 7r..7r] and oscillate around 6' = ( Fig. [^Jj). In this 
regime the population oscillation amplitude is expected 
to decrease with decreasing q, i.e. at lower magnetic field, 
as less population transfer is necessary to create a mean 
field energy of the size of the quadratic Zeeman shift. 
This decrease in amplitude is also confirmed by the data 
in Fig.Hl 

In the resonance region at intermediate magnetic field, 
the relative movement of the (0) and (±1) waves is very 
slow, with 9 being on the border of oscillation and con- 
tinuous increase. This leads to long time unidirectional 
spin mixing and thus maximum amplitude. 

We want to point out that also the initial direction 
of spin mixing increasing the (0) wave amplitude as ob- 
served in Fig. ^ can be explained by simple arguments. 
The chosen initial state C'"' has zero phase shift = 
and shows no mean field dephasing between the (±1) and 
(0) waves. This is due to the fact, that it originates from 
the fully stretched (—2) eigenstate subjected to a 90 de- 
gree rotation). As the mean field part of the Hamiltonian 
is invariant under rotations C'"' stays a mean field eigen- 
state. Consequently the quadratic Zeeman energy shift 
will always cause the (±1) wave to evolve faster then the 
(0) wave or in other words 9 will initially always grow. 
The transfer of particles going in the direction of the wave 
lagging behind will thus result in an initial increase of the 
(0) wave (see Fig.|21). 

In addition to the understanding of the tendencies ob- 
served in Fig. ^ within the single mode approximation, 
a quantitative comparison allows to deduce the validity 
limits of this approximation for the elongated trap ge- 
ometry used in the experiment. The breakdown of the 
single mode approximation is accompanied by the strong 
damping of the oscillation at low q, which is probably due 
to the onset of spatial structure formation in spin dynam- 
ics jlOj . In contrast to the quadratic Zeeman dominated 
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FIG. 3: TOP: Deviation (rms of ttif = population) of 
SMA theory fit {gi{n) = 54 s~-^ fixed) from our data (Fig.^J. 
Tfie abscissa specifies tfie time range over which the fit was 
perfomed. BOTTOM: two examples of fits compared to data. 



regime spin dynamics in tlie mean field regime is charac- 
terized by a density dependent coupling constant gi{n). 
For mean field driven spin dynamics the higher density 
parts of the spinor condensate will thus dephase relative 



to the lower density parts and population oscillations in 
the total fractions will be washed out. 

The validity range of the SMA can also be inferred from 
our experimental data (Fig. |2J). While the evolution of 
Spinor condensates in our elongated geometry starts to 
deviate from SMA behavior after « 5 ms in the mean 
field regime, there is a near perfect agreement over more 
than 20 ms in the mean field region. Our most important 
observation in this respect is that even for condensate 
sizes much above the spin healing length the single mode 
approximation works remarkably well once the quadratic 
Zeeman effect dominates the relative phase evolution. 

In conclusion we have investigated a magnetically tun- 
able resonance in coherent spin dynamics and have ana- 
lyzed this phenomenon in terms of phase matching con- 
siderations, which are applicable to general values of to- 
tal spin. In the future this resonance will allow the con- 
trolled use of spin dynamics to fully convert one spinor 
state into another, even with local control by shifting in 
and out of the validity range of the single mode approx- 
imation. These possibilities open new perspectives for 
spatiotemporal tailoring of the spinor condensate wave- 
function, e.g. for the creation of highly nonclassical spin 
states with spatially varying degree of entanglement. 
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